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Abstract 

Initial state radiation (ISR) in deep inelastic scattering of charged leptons 
on nucleons can lead to substantial shifts in the reconstruction of and the 
scaling variables x and y. By measuring at least three out of the four variables: 
energy and production angle of the lepton and the hadron system, the energy 
of the ISR photon emitted under zero degrees can be determined. The ISR 
effect is illustrated for experimentation at HERA with different methods of 
kinematic reconstruction. 

1 Introduction 

The kinematics of deep-inelastic electron proton scattering (DIS) by neutral current 
exchange (NC), 

ep^eX, (1) 



is characterized by the mass squared (— Q^) of the exchanged current and the Bjorken 
scaling variable x. In a hermetic detector, and x can be determined from the 



energy and angle of the scattered electron or from the energy and momentum com- 
ponents of the produced hadron system or from a combination of any other two of 
these quantities 0, 

The determination of x and is erroneous when the incoming electron radiates a 
photon before the hard interaction takes place (initial state radiation, ISR) Q. When 
averaged over a large DIS event sample, the effect of ISR in general is rather small 
0. However, ISR can have an important effect when, due to ISR, events produced in 
a — X region of large cross section migrate into a region where the cross section is 
expected to be small. An example for such a constellation is presented by the possible 
excess of high-x, high-Q^ events recently observed at HERA ^f. 

Most of the ISR photons are emitted very close to the direction of the electron such 
that the emission angle can be assumed to be zero (peaking approximation). The ISR 
energy and the Q^, x values can then be obtained on an event-by-event basis from a 
measurement of three out of the four quantities mentioned above ^, ^, Q or from 
a global fit [0]. 

In charged current (CC) processes, 

ep-^uX, (2) 

where the u escapes detection the information from the hadron system can be used 
to put a limit on the energy emitted by ISR while a measurement requires direct 
observation of the ISR photon. 

This note presents the relevant formulae in a form suitable for HERA experiments 
and illustrates the migration of events due to ISR in the x — y,x — Q'^,y — Q'^ and 
M — y planes, where M is the effective mass of the electron current quark system. 



2 DIS kinematics in the absence of ISR 

In deep-inelastic ep scattering (eq. |l]) depicted in Fig. ^ the incoming electron couples 
to the neutral electroweak current J = 'j, Z^. The relevant variables in the absence 
of QED radiative effects are: 



^Due to the comparatively large mass of the proton, ISR from the incoming proton can be 
neglected. In principle, photon emission by the scattered electron (final state radiation, FSR) will 
affect also the determination of the kinematics; in practice, if the scattered electron is detected in a 
calorimeter, the effect will be small because the signals from the electron and the FSR photon are 
summed. 

^For DIS at HERA with between 10 and 100 GeV^, roughly 6% of the events have an ISR 
photon carrying between 5 and 20% of the energy of the incident electron. 
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Ep proton beam energy 

Eg electron beam energy 

p = (0, 0, Ep, Ep) four momentum of incoming proton with mass rrip 

e = (0, 0, —Eg, Eg) four momentum of incoming electron 

e' = (EgSinO'g, 0, E^cosOg, E'^) four momentum of scattered electron 

s = (e + pY = 4:EpEf, square of total ep cm. energy 

= (e — e'Y = —Q^ mass squared of exchanged current J 

= square of four momentum transfer 

V = q - p/ iTLp energy transfer by J in p rest system 

T-'max = s/{2mp) maximum energy transfer 

y = {q ■ p)/{e ■ p) = v j Vmax fraction of energy transfer 

X = Q"^/ (2g ■ p) = j (ys) Bjorken scaling variable 

qc = X ■ p + {e — e') four momentum of current quark 

= {e' + qc)"^ = X ■ s mass squared of electron - current quark system. 



The electron and proton masses {me,mp) have been neglected. The angle 6'^ of the 
scattered electron is measured with respect to the incoming proton (see Fig. |l]b). For 
simplicity, the scattered electron has been assumed to have zero azimuthal angle. For 
fixed cm. energy a/s, the inclusive process ep — >• eX is described by two Lorentz- 
invariant variables which can be determined from the electron or the hadronic system. 

The electron side yields: 

y = I- ^(1 - cose'g) 

Q2 = 2E,E'^{1 + cosO'J 

X = E'g{l + cose',)/{2yEp) 
E'g = {I - y)Ee + xyEp 

cose'g = ^yEp-i^-y)Ee 

xyEp + (1 - y)Ee 
E'/sinX = Axy{l-y)EgEp. (3) 



With the help of the relation M = y/x ■ s the variable x can be replaced by the mass 
M of the electron current quark system. 

Neglecting the proton remnant, the hadron system j, with energy Ej and production 
angle 9j, yields: 



E^ 
2E, 

12 Z7'2„-„2 



Q' = E'sin%/il - y) 



Ej = yE, + x{l - y)Ep 

-yEe + {1- y)xEp 

cosOj = — ^ — — ^ 

^ yEe + {l-y)xEp 

E]stn^e^ = Axy{l-y)E,Ep = Q\l-y). (4) 

Using the method of Jacquet-Blondel |^ the hadron variables can be determined 
approximately by summing the energies {Eh) and transverse {pTh) and longitudinal 
momenta {pzh) of all final state particles. The method rests on the assumption 
that the total transverse momentum carried by those hadrons which escape detection 
through the beam hole in the proton direction as well as the energy carried by particles 
escaping through the beam hole in the electron direction can be neglected. The result 
is: 

yjB = 7^J2(^h~ PZh) 



(EhPxh)^ + (EhPYh) 
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Q^jB — 1 

1 - yjB 

xjB = Q^bUvjbs) (5) 



The double-angle (DA) method of ^ uses the electron scattering angle and the 
angle 7/1 which characterizes the longitudinal and transverse momentum flow of the 
hadronic system (in the naive quark-parton model 7/1 is the scattering angle of the 
struck quark): 

^ {HhPXh? + {HhPYhf - {HhjEh - PZh)? ( . 



leading to 



g2 _ AElsin-ih (1 + cosO'^) 

Ee siri'jh + sinO'^ + sin{'yh + 0'^ 

Ep siri'^h + sinO'^ — sin{'^h + O'e 
siwyh (1 + cos6'^ 



xda = 
Vda = 



_ (1 - cos'^h) sinO'^ 

sifi'jh + sin9'^ — sin{pfh + Q'e) 
= QIaUxdas). (7) 



4 



The energy of the scattered electron obtained by the DA-method is: 



' siwjh + sin9'^ — sini^jh + ^e) 



The S method [0| combines the measurements of the energies and angles of the 
electron and hadron system. Making use of the relation 

J2iEh-pzh) + E'^{1 - cosO'J = 2Ee, (9) 

h 

which holds exactly for a hermetic detector, one can write 

J2hiEh-pzh) 



2/s 

Ql 



Eh{Eh-pzh) + E>{l - cosO',: 
E' sin^e' 



1 - 2/s 



X. = ^. (10) 



3 DIS kinematics with ISR 

The photon from ISR is assumed to be emitted in the direction of the beam electron 
(see Fig. ^). In order to avoid confusion, the kinematic quantities of the beam particles 
will be denoted by the subscript beam] those of the particles participating in the 
hard interaction will be denoted by the subscript true. The "measured" variables 
are those obtained under the assumption of no ISR emission. They will be identified 
by the method of reconstruction (i.e. e, DA or JB method) and by the subscript 
meas. They will be assumed to be measured without errors. For ease of reference all 
relevant variables are listed below: 



Initial state (nominal): 



Epbeam i -^ebeam 

^nominal ^Epiygd^n E^ij^am 



proton and electron beam energies 
square of nominal ep cm. energy 



ISR photon: 



E^ 

Eetrue (1 /'7)-^ebeam 



energy of the ISR photon 

fraction of electron beam energy carried by ISR 
energy of incoming electron in the hard interaction 
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Hard interaction: 



•^true 1 Vtrue i Q 



emeasy Vemeasy ^ emeas 
e emeas ' emeas 
•^DAmeasi VDAmeasi Q DAmeas 

w 

e DAmeas 
^JBmeasi lljBmeasi QjBmeas 



true 



xe, ?/s, Ql 



true values 

X, y, obtained with the e-method 
E'^, 6'g obtained with the e-method 
X, y, obtained with the DA method 
E'^ obtained with the DA-method 
X, y, obtained with the JB method 
X, y, obtained with the S method 



The effect of ISR is taken into account by replacing in eqs. |^ to the electron beam 
energy by E^true- The result is 



•Sir 



nominal 



'IV. 



and for the electron method: 



Xfr 



Vtr 



Q 



true 



Ve 



Vtrue 
y emeas (1 /' 



7/ 



^emeas 



y emeas 

- A 



1 - A 

2 

emeas 



Q 



1 - A) • 



'121 



One finds for the electron method that 



cannot exceed y emeas- A < 
Xtrue < 1 which Icads to < ^ 



7 ^ y emeas- J 7 

y emeas (1 -^emegs) 



. f-y is limited also by the condition 



emeas y emeas 



all three variables (x, y, Q^) are affected by ISR. In the presence of ISR the 
points move in the x, y plane along lines of constant Qlmeas- note, XemeasVemeas = 



)2 

' emeas 



XtrueVtrue and therefore Q 

ISR reduces s the true value is smaller than Q 



^nominal emeas V emeas ^nominal -^true Vtrue - 
2 

emeas ' 



Since 



the measured x value is smaller than the true one, x emeas < Xtrue] 



the measured y value is larger than the true one, y emeas > Vtrue- At small y the 
shift in y - and therefore also in x - can become quite large, e.g. for L = 0.1, 



Vtr 



0.2 one finds ye 



I. Ay true and Xe 



0.7x 



true- 
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The corresponding results for the DA variables are: 

•^true \^ J-y) ' •^DAmeas 

(13) 

Vtrue y D Ameas 

Qtrue ~ QoAmeas ' (-'- ~ fl) ■ 

One finds for the DA method that 

• X and are affected by ISR while y is insensitive to ISR; 

• the measured x value is larger than the true one, xd Ameas > Xtme] 

• ISR has a strong effect on Q"^ since the factor (1 — /^) enters quadratically. To 
give an example, for /-y = 0.2 and Ql^ue = 20000 GeV^, the measured value 
is more than 50% larger, viz. Qfj^^eas = 31250 GeV^. 

The results for the Jacquet-Blondel method are: 

1 y J Bmeas ( ^ A\ 

•^true ■^JBmeasZ 77Z r \ \^^) 

J- — y J Bmeas /[^ " Jj) 

1 

ytrue y J Bmeas Z 
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^true ^ J Bmeas i 



y J Bmeas /{^ — f-y) 

One finds for the JB method that 

• cannot exceed (1 - yj Bmeas)- f'y<i)- - yj Bmeas)- fy is limited also by the 
condition Xtme < 1 which leads to L < (]-''JBrueas)(i-yjB^.as) 

• x,y and are affected by ISR. 
The results for the S method are: 

_ XT,meas c-x 

•^true — r \^^) 

ytrue = 2/S 
Qtrue ~ Qt,- 

One finds for the S method that 

• y and are not affected by ISR. This is so because eq. ^ holds also in the 
presence of ISR when E^true is used for E^, instead of E^beam- 

• Note that M is also unaffected by ISR since = Xtrue^true = Xj^measS nominal - 
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3.1 Graphical illustration of ISR effects 



The effects of ISR on the measurements of x, y and when using the e, DA, JB 
and S methods have been calculated for colhsions of Eebeam = 27.5 GeV electrons on 
Epbeam = 820 GcV protous. The results for the x — y,x — Q'^,y — Q"^ and M — y 
planes are presented in Figs. |^ to |16| for two values of the fractional energy carried 
away by the radiated photon: = 0.1 and 0.2. The arrows shown point from the 
measured to the true quantities. The curves in the x — y and M — y diagrams show 
lines with constant Ql 



)2 

/ meas • 



3.2 Determination of the energy of the ISR photon 

In selecting DIS events the HERA experiments HI and ZEUS place a cut on the 
difference between the total energy and longitudinal momentum 

6 = Y.{E- Pz) (16) 

observed in the detector which puts a limit on the energy of possible ISR. For a 
hermetic detector 6 = 2Eebeam- If the ISR photon escapes detection the requirement 
5 > ^min limits the photon energy: 

and therefore 

fj < ^ — Smin/i'^Eheam), (18) 

provided that measurement errors and other detector effects can be neglected. 

ISR in general affects the variables x, y, differently when measured with the elec- 
tron, DA or JB methods. For instance, in the presence of ISR the x value measured 
with the electron method is always smaller than the true value, Xtme, while the x 
value obtained with the DA method is always larger than Xtme- Turning the ar- 
gument around, if one observes Xemeas < xuAmeas and Qlrneas < QnAmeas then this 
is a strong indication for the emission of ISR, provided measurement errors can be 
neglected. 

The energy carried by the ISR photon can be determined directly from a comparison 
of the variables found with the different methods. Assuming that measurement errors 
can be neglected one obtains 



8 



from the comparison of e and DA variables: 



E' 

£ -l eemeas f -\ r\\ 

h - ^ - y, 

e DAmeas 

I -, •^emeas \ 

Vemeas \^ ) 
■^DAmeas 

Vemeas UDAmeas 



1 - Vda 

Q: 



meas 

2 

emeas 



QoAm eas 

from the comparison of e and JB variables: 

f-y Vemeas VJBmeas (^0) 

f-i „, '^JBmeas\ 
— VJBmeas )[^ 



emeas 



and from the comparison of DA and JB variables: 



J ^ VjBmeas (21) 

V D Ameas 

/-, \ /T -VjBmeas \ 

= (1 - VJBmeasjyi^ ) 

•^DAmeas 

In the S method the result is 

Of course, the most precise result for NC processes on can be expected from a 
common fit to all variables measured as advocated in ||^. 

As mentioned before, for CC processes, the ISR energy must be determined from the 
direct observation of the photon. The information from the hadron side allows only 
to put a limit on the photon energy: < (1 — VJBmeas)- 
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Figure 1: DIS, ep — > eX, without ISR: diagram (top) and event configuration in a 
HERA detector (bottom). 
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^current 




Figure 2: DIS, ep — > eX, with ISR: diagram (top) and event configuration in a HERA 
detector (bottom). 
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X 



Figure 3: The x — y plane: The effect of ISR on the determination of x and y by 
the electron method for = — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {xemeas,yemeas) rcconstructcd with the electron method to the point 
{x true, y true) ■ The curves show lines with constant Ql^neas- 
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Figure 4: The x — y plane: The effect of ISR on the determination of x and y by 
the DA method for L = -^^^ — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {xDAmeasiU DAmeas) reconstructed with the DA method to the point 
{x true, y true) ■ The curves show Unes with constant QoAmeas- 
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Figure 5: The x — y plane: The effect of ISR on the determination of x and y by the 
JB method for L = -e;^ — = 0.1 (top) and 0.2 (bottom). The arrows point from the 
point {x J Bmeas,y JBmeas) reconstructed with the JB method to the point {xtme, Vtrue) ■ 
The curves show lines with constant Q^jBmeas- 



15 



=0.1 (x,y)5,^,,, ^ (x,y) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



=0.2 (x,y),^^^^ ^ (x,y), 




0^ 

^ Tmeas 

.bOGOC'GeV" 

.40G00GeV' 

.-30000GeV' 
:25000GeV' 
-20000GeV' 
SOOOGeV' 

----.-^ 1 OOOOGeV' 
SOOOGeV' 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 6: The x — y plane: The effect of ISR on the determination of x and y by 
the S method for L = t?^^' =0.1 (top) and 0.2 (bottom). The arrows point from 
the point (aj^meas! I/Emeas) reconstructed with the E method to the point {xtme, Vtrue)- 
The curves show lines with constant Qj^^^as- 
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Figure 7: The x — plane: The effect of ISR on the determination of x and by 
the electron method for = — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {xemeas, Qlmeas) reconstructed with the electron method to the point 

rue-, Qtrue)- 
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Figure 8: The x — plane: The effect of ISR on the determination of x and 
by the DA method for = — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {xDAmeas,QDAmeas) reconstructed with the DA method to the point 
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Figure 9: The x — plane: The effect of ISR on the determination of x and 
by the JB method for L = — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {xjBmeas,QjBmeas) reconstructed with the JB method to the point 

rue, Qtrue)- 
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Figure 10: The x — Q"^ plane: The effect of ISR on the determination of x and 
by the S method for = — = 0.1 (top) and 0.2 (bottom). The arrows 
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Figure 11: The y — Q"^ plane: The effect of ISR on the determination of y and by 
the electron method for = — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {y emeas, Qlmeas) reconstructed with the electron method to the point 
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Figure 12: The y — plane: The effect of ISR on the determination of y and 
by the DA method for = — = 0.1 (top) and 0.2 (bottom). The arrows point 
from the point {yoAmeas, QoAmeas) reconstructed with the DA method to the point 

rue, Qtrue)- 
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Figure 13: The y — plane: The effect of ISR on the determination of y and 
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Figure 14: The M — y plane: The effect of ISR on the determination of M = y/xs 
and y by the electron method for = — = 0.1 (top) and 0.2 (bottom). The 
arrows point from the point {Mgrneas,y emeas) reconstructed with the electron method 
to the point {Mtrue^Vtrue)- The curves show lines with constant Qlrneas- 
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Figure 15: The M — y plane: The effect of ISR on the determination of M = y/xs 
and y by the DA method for L = -et^ — = 0.1 (top) and 0.2 (bottom). The arrows 
point from the point {MDAmeas,y DAmeas) reconstructed with the DA method to the 
point {Mtrue, Utrue) ■ The curves show lines with constant QDAmeas- 



25 



f =0.1 (M,y)„„^, ^ (M,y)^ 




eV 
eV 



50 75 100 125 150 175 200 225 250 275 300 



^1 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




M(GeV) 
=0.2 (M,y)jB^^^^ ^ (M,y)^^^ 




I I I I I I I I 



^ JBmeos = 

SOOOOGeV' 

40000GeV" 

OOOGeV' 
OOOGeV" 
OOOGeV' 
OOOGeV' 
COOGeV' 

mm' 



50 75 100125150175 200225 250275 300 



M(GeV) 



Figure 16: The M — y plane: The effect of ISR on the determination of M = y/xs 
and y by the JB method for L = -^^^ — = 0.1 (top) and 0.2 (bottom). The arrows 
point from the point {M j smeas ^ U j Bmeas) reconstructed with the JB method to the 
point {Mtrue, Utrue) ■ The curves show lines with constant Q^jBmeas- 
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